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We present the Standard Model prediction for the eight angular coefficients of the W boson, which 
completely describe its differential cross section in hadron collisions. These coefficients are ratios of 
the W helicity cross sections and the total unpolarized cross section. We also suggest a technique 
to experimentally extract the coefficients. 


I. INTRODUCTION 


The analytical study of the W boson is essential for 
the understanding of many open questions related to the 
electroweak physics, like the origin of the electroweak 
symmetry breaking and the source of the CP violation. 
Since its discovery the W hadronic cross section, mass, 
and width have been measured with great precision [ 1 ], 
On the other hand complete physical information is con¬ 
tained in the boson’s angular distribution in three di¬ 
mensions, given by its differential cross section, which 
can be written as a sum of helicity cross sections. These 
quantities are related to the nature of the electroweak 
processes, the W polarization, and the presence of QCD 
effects. In this paper we address these issues in the case 
of W produced in hadron collisions. 

The total differential cross section of the W production 
in a hadron collider [ 2 ] is given by the equation: 

da _ 3 da u + ^ 

dq^dyd cos 9d<f> 167 r dq^dy 

+ -Ao(l — 3 cos 2 9) + Ai sin 20 cos</> 

+ ^A 2 sin 2 9cos 2(f> T A 3 sin# cos ^ 

+ A 4 cos 9 + A 5 sin 2 9 sin 2 (f> 

+ A e sin 29 sin (f> + Ay sin 9 sin (j)\ (1) 


where qr and y are the transverse momentum and the ra¬ 
pidity of the W in the lab frame and 9 and </> are the polar 
and azimuthal angles of the charged lepton from the W 
decay in the Collins-Soper (CS) frame [3]. The CS frame 
is used because in this frame we can experimentally re¬ 
construct the azimuthal angle </> and the polar quantity 
| cos0|. Our ignorance of the W longitudinal momentum, 
which is due to our inability to measure the longitudinal 
momentum of the neutrino in a hadron collider, intro¬ 
duces a two-fold ambiguity on the sign of cos 9. The 
quantity da u /dq^dy is the angles-integrated unpolarized 
cross section. 

The dependence of the cross section on the leptonic 
variables 9 and (f> is completely manifest and the depen¬ 
dence on the hadronic variables cpf and y is completely 


hidden in the angular coefficients Ai(qT,y). This allows 
us to treat the problem in a model independent man¬ 
ner since all the hadronic physics is described implic¬ 
itly by the angular coefficients and it is decoupled from 
the well understood leptonic physics. The angular coeffi¬ 
cients are ratios of the helicity cross sections of the W and 
da u /dq^dy. In order to explicitly separate the hadronic 
from the leptonic variables, the helicity amplitudes were 
used to describe the hadronic tensor associated with the 
hadronic production of the W [4]. The leptonic tensor on 
the other hand is analytically known, leading to analytic 
functions of the angles of the charged lepton in Equa¬ 
tion (1). It is common to integrate Equation (1) over y 
and study the variation of the angular coefficients as a 
function of qr- 

If the W is produced with no transverse momentum, 
it is polarized along the beam axis because of the V-A 
nature of the weak interactions and helicity conservation. 
In that case A 4 is the only non-zero coefficient. If only 
valence quarks contributed to the W T production, A 4 
would equal 2 , and the angular distribution ( 1 ) would be 
~ (1 ± cos 9) 2 , a result that was first verified by the UA1 
experiment [5]. 

If the W is produced with non-negligible transverse 
momentum, balanced by the associated production of 
jets, the rest of the angular coefficients are present and 
the cross section depends on the azimuthal angle 4 > as 
well. The last three angular coefficients - A 5 , Aq , and 
A 7 - are non-zero only if gluon loops are present in the 
production of the W [ 6 ]. Hence, in order to study all 
the angular coefficients and associated helicity cross sec¬ 
tions of the IT in a hadron collider, we have to consider 
the production of the W with QCD effects at least up to 
order a 2 . 

The importance of the determination of the W angu¬ 
lar coefficients is discussed in [7], and is summarized here. 
This study allows us to measure for the first time the dif¬ 
ferential cross section of the W and study its polarization, 
since the angular coefficients are related to the helicity 
cross sections. It also helps us verify the QCD effects in 
the production of the W up to order a 2 . In addition, 
A 3 is only affected by the gluon-quark interaction and 
its measurement could constrain the gluon parton dis- 
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FIG. 1: The charged lepton <fi distribution in the Collins- 
Soper W rest-frame for four qr regions. The distributions are 
normalized to unity. 

tribution functions. Moreover, the next-to-leading order 
(NLO) coefficients A 5 , A 6 , and A 7 are P-odd and T-odd 
and may play an important role in direct CP violation 
effects in W production and decay [ 6 ]. Finally, a quanti¬ 
tative understanding of the W angular distribution could 
be used to test new theoretical models and to facilitate 
new discoveries. 

In this paper we present the Standard Model predic¬ 
tion for the angular coefficients A* as a function of qr, 
for proton-antiproton collisions at 1.8 TeV, using the 
DYRAD Monte Carlo program [ 8 ], an event generator 
of W+ jet up to order a 2 . Three of the coefficients, in¬ 
cluding A 4 which is the dominant one up to qt = 100 
GeV, are presented for the first time. We also suggest a 
new method for the extraction of the angular coefficients. 


II. STANDARD MODEL PREDICTION FOR 
THE W ANGULAR COEFFICIENTS 

In order to study the angular distribution of the W we 
have to choose a particular charge for the boson. In this 
paper we present the results for the W~. The angular 
coefficients for the W + can be extracted by CP trans¬ 
formation. In the Collins-Soper frame, the CP transfor¬ 
mation leaves </> unchanged and takes 9 to 7 r — 9. If we 
assume that Equation (1) describes W~ bosons, we have 
to change the sign of coefficients A\. A 4 , and A§, in order 
to describe W + bosons, without changing the definition 
of the Collins-Soper frame. 

We generate Monte Carlo IF+jet events up to a 2 , in¬ 
cluding up to one gluon loop, from proton-antiproton 
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FIG. 2: The charged lepton cos# distribution in the Collins- 
Soper W rest-frame for four qr regions. The distributions are 
normalized to unity. 


collisions at y/s = 1.8 TeV, using the DYRAD gen¬ 
erator. We run with minimal kinematic and accep¬ 
tance cuts, with a minimum transverse energy for the 
je t of 10 GeV and jet -jet angular separation of A R = 
-\/(A? 7 i a b ) 2 + (A^i a b ) 2 = 0.7, in the pseudorapidity-phi 
space in the lab frame. The CTEQ4M parton distribu¬ 
tion functions were used. 

We measure the 9 and 0 angles of the charged lep¬ 
ton in the Collins-Soper W rest-frame. At the Monte 
Carlo event-generator level, we know the momentum of 
the neutrino, so there is no two-fold ambiguity on the sign 
of cos 9. The CS frame is the rest-frame of the W where 
the z-axis bisects the angle between the proton momen¬ 
tum (pcs) and the opposite of the antiproton momentum 
(—pcs) in the CS frame. The signs of the angular coef¬ 
ficients depend on the way the Collins-Soper a;-axis and 
p-axis are defined. In this paper, we define them so that 
the x-z plane coincides with the pcs — PCS plane and the 
positive p-axis has the same direction as pcs x PCS- The 
SM distribution of the 0 and cos 9 for four qr bins (15-25, 
25-35, 35-65 and 65-105 GeV) is shown in Figures 1 and 
2 respectively. We note that at low transverse momen¬ 
tum of the W, the 0 distribution is almost flat, whereas 
the cos# distribution almost follows the (1 + cos Q ) 2 law. 
In the fourth qr bin, there is a strong ^-dependence of 
the cross section and the cos 9 distribution is almost a 
straight line (|cos#| is flat). There is a correlation be¬ 
tween cos 9 and A , with low cos 9 corresponding to low 
events. 

To calculate the angular coefficients from the angles of 
the charged lepton, we use the method of moments. We 
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FIG. 3: The Standard Model prediction for the angular coef¬ 
ficients of the IF produced in a collider at y / s=1.8 TeV. QCD 
effects are included up to order a 2 s . 


first define the moment of a function m( 6 , <j>) as 

/ / da(q T , y , 9, <t>)m(9, (j>)d cos 6 d(j> 
f f da(qx,y, 9, <f>)dcos 9d<j> 

We can easily prove that: 


( 2 ) 


< too > 

< m i > 

< m .2 > 

< m3 > 

< m 4 > 

< m3 > 

< m.Q > 

< m 7 > 


<l(l-3cos 2 0) >= 

< sin 29 cos (j) >= —A\ 

5 

< sin 2 9 cos 2(j> >= — A 2 

< sin 9 cos <f> >= - A 3 

< cos 9 >= -A 4 

< sin 2 9 sin 2 / >= —A 5 

5 

< sin 29 sin <p >= -A e 

5 

< sin 9 sin <f> >= —Ay 


( 3 ) 

For a set of discrete generator (or experimental) data, 
we substitute the integrals of Equation (2) by sums and 
the cross section values by the weights w t of the Monte 
Carlo events. 


< m(9, (j)) >= 


J2iLi m (6iAi)wi 

Ei=1 W i 


( 4 ) 


FIG. 4: The Standard Model prediction for the next-to- 
leading order coefficients (As, Aq, and Ay) and Ai for the 
W production in a collider at \/s=1.8 TeV. The coefficients 
A 5 , Aq, and Ay are present only if gluon loops are included 
in the calculation and they are T-odd and P-odd. 


By solving Equations (3) for the angular coefficients 
and substituting the moments by the discrete expres¬ 
sions (4), we extract the Standard Model prediction. By 
ignoring the W rapidity , we actually calculate the y- 
integrated angular coefficients, which are now functions 
of just q t- The results are shown in Figure 3. The an¬ 
gular coefficients Ay, A 4 , and A 6 are presented for the 
first time. The DYR.AD Monte Carlo generator is more 
reliable for qt > 10 GeV, which is also the transverse 
energy cut for our jets, and this value determines the 
minimum of our qx- axis. The maximum is determined 
by the Monte Carlo statistics. 

We notice that indeed A 4 is the only surviving major 
coefficient at low qt values. It is also the only angular co¬ 
efficient that decreases as qr increases. The angular coef¬ 
ficient Ai, although it is a leading order (LO) coefficient, 
is much smaller than the leading order coefficients (Ao, 
A 2, A 3 , and A 4 ) and comparable to the next-to-leading 
order ones (A 3 , Aq, and Ay). The coefficients Aq and A 2 
would be exactly equal if gluon loops were not included 
[2], At order a 2 , A 0 is consistently greater than A 2 . To 
better determine A±, A 3 , Aq, and Ay we use the four qx 
bins, to improve the statistics. The result is shown in 
Figure 4. There are relations that directly connect the 
angular coefficients with the helicity cross sections of the 
W [2]. We first extract the unpolarized cross section of 
the IF as a function of qx and using our prediction for 
the angular coefficients, we arrive at the Standard Model 
prediction for the IF helicity cross sections at 3 /s = 1.8 
TeV, shown in Figure 5. Here dcq is the helicity cross 
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FIG. 5: The Standard Model prediction for the helicity cross 
sections of the W produced in a collider at ^/s=1.8 TeV. QCD 
effects are included up to order a 2 a . 


FIG. 6: The experimentally observed electron <j> distribution 
in the Collins-Soper W rest-frame for four qr regions. The 
distributions are normalized to unity. 


section that corresponds to the angular coefficient 


III. EXPERIMENTAL DETERMINATION OF 
THE W ANGULAR DISTRIBUTION 


E. Mirkes [2] first realized the problem of directly mea¬ 
suring the angular coefficients. The angular distributions 
of Figures 1 and 2 are seriously distorted after the ef¬ 
fects of the detector are considered and quality cuts are 
imposed on the data sample. To study the effect, we 
treat the generator leptons as electrons and we pass them 
through a detector simulator [9]. The new </> and cos 9 
distributions are shown in Figures 6 and 7 respectively. 
The shapes of the muon distributions are the same as 
that for electron W+ jet events, however less events are 
detected, because of the lower muon acceptance of a typ¬ 
ical hadron collider detector. The main reason for the 
difference between Figures 1, 2 and Figures 6, 7 is the 
leptons transverse momentum cuts (p T > 20 GeV and 
Pj, > 20 GeV) and the charged lepton rapidity cut (cen¬ 
tral leptons are considered, \y l \ < 1). 

The problem of distortion of the (f> and cos# distri¬ 
butions due the detector effects and quality cuts is not 
the only one. A more fundamental problem is the actual 
measurement of these angles. To measure them, we need 
to reconstruct the W in the three dimensional momen¬ 
tum space, in order to boost to its center of mass. The 
longitudinal momentum of the neutrino is not measured, 
but it is constrained by the mass of the W. based on 






Electron cos6 (35 < p T w < 65 GeV) Electron cos0 (65 < p T w < 105 GeV) 


FIG. 7: The experimentally observed muon cos# distribution 
in the Collins-Soper W rest-frame for four qr regions. The 
distributions are normalized to unity. 


equation: 

Pz = y [a P 1 z ± E l \JA 2 — 4(py) 2 (py) 2 ^ (5) 

where 

A = M^y + q? — {p l T ) 2 — {pt ) 2 ) 












IV. EXPERIMENTAL EXTRACTION OF THE 
ANGULAR COEFFICIENTS 
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FIG. 8: The error in <j> reconstruction due to the selection 
of the W mass is considerably lower than the error in cos 9 
reconstruction, especially at high Mx' values (upper plots). 
This corresponds to better reconstruction of the <j>, while the 
cos 9 reconstruction is worse, especially at low values of | cos 9 \ 
(lower plots). 


E l is the energy of the charged lepton, p l T and are 
the transverse momentum of the charged lepton and the 
neutrino and p l z is the longitudinal momentum of the 
charged lepton. The two solutions for the longitudinal 
momentum of the neutrino lead to two solutions for the 
W longitudinal momentum. Both solutions correspond 
to the same <f> but to opposite cos 8 values. 


Moreover, according to Equation (5), for each event, 
we have to input a mass for the W to get p v z and even¬ 
tually <j> and cos 6 \. The mass of W is not known on 
event-by-event basis, we just know its pole mass and its 
Breit-Wigner width. Based on these two established val¬ 
ues, we can plot the uncertainty in the measurement of 
<f> and | cos0| introduced by the uncertainty in the mass 
of the W. For each Monte Carlo event we generate W 
masses that are greater than the transverse mass for the 
particular event and follow the Breit-Wigner distribu¬ 
tion. In Figure 8 we see that the systematic error on the 
measurement of <f> is very small, but the cos 9 system¬ 
atic error is significant, especially at low | cos 01 and at 
big values of the transverse mass of the W. This makes 
the direct measurement of the | cos 0| distribution more 
challenging. 


In [7] it is suggested that the experimental distribu¬ 
tions of Figures 6 and 7 should be divided by the Monte 
Carlo distributions obtained using isotropic W decays. 
This method results in distributions similar to those 
shown in Figures 1 and 2 and the extraction of the angu¬ 
lar coefficients is easier. Here we present a method that 
does not bias the experimental data by Monte Carlo data. 
Instead, it uses the knowledge of the detector and its ef¬ 
fect on the theoretical distributions. We will demonstrate 
the method for the <p analysis. 

If we integrate Equation (1) over cosf? and y, we get: 

, . = C"(l + Pi cos <j> + p 2 cos 2 (j> 

dq^dtp 

+/?3 sin <f> + Pi sin 2<j>) (6) 


where 


r > _ 1 da R _ 37r a R _ a 2 
„ 3-7T Ak 

P 3 = A r , Pa = -y 


(7) 


The observed <j> distribution is given by Equation (6), 
only if we ignore the effects of the detector and kine¬ 
matic cuts. In any other case, there is an acceptance 
and efficiency function ae(qr, cos 0, </>) which multiplies 
(1) before it is integrated over cos0 and as a result, no 
angular coefficient is completely integrated out. In the 
actual data, what we measure is the number of events, 
which is: 


N(q T , <t>) 


/ d(j 

-—— --ae(g T , cos 9, p)d cos 9 

dqrdcpd cos 9 

+ Nb g (q T , 4>) 


Cdt 

( 8 ) 


where C is the luminosity and ae(qr, cos 9, (f>) is the ac¬ 
ceptances and efficiencies for the particular W trans¬ 
verse momentum and pixel in the (cos 9, (f>) phase space. 
Nbgiqr , p) is the background for the given <f> bin and qr- 
If we combine Equations (8) and (1), then the measured 
distribution is: 


7 

N(q T , P) = C'{f -1 + A-ifi) + N bg (q T , P) (9) 
1=0 


where C' = C f Cdt , and /) are the fitting functions, 
integrals of the product of the explicit functions of cos 9 
and p and the ae(cos0, p): 


PIT 

fi{q T ,<j>)= / (j>)ae(q T , cos0, c/))dcos 0, (10) 

Jo 

.,7 


i = — 1, 
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FIG. 9: The leading order f-i,--. , 4 ( 0 ). These functions are 
multiplied by the respective W angular coefficients to give us 
the experimentally observable 0 distributions. 


FIG. 10: The next-to-leading order/5 : 6,7(0). These functions 
are multiplied by the respective W angular coefficients to give 
us the next-to-leading order corrections to experimentally ob¬ 
servable 0 distributions. 


where 

9-i(9,0) 

gi(9,0) 

93 ( 9 , 0 ) 

95 ( 9 , 0 ) 

97(9,0) 


1 + cos 2 9, 
sin 20 cos</>, 
sin 9 cos 0 , 
sin 2 9 sin 20, 
sin 9 sin 0 


go(9, 0) = 1/2(1 - 3 cos 2 9) 
92 (9,0) = 1/2 sin 2 0 cos 20 
54 ( 0 , 0) = cos 9 
ge(9,0) = sin 29 sin 0 


The fi functions can be calculated explicitly if we know 
the acceptance and the efficiency of the detector. Be¬ 
cause we multiply by ae(qr, cos 9, 0) before integrating 
over cos0, no fi is exactly zero. As a result, all coeffi¬ 
cients are in principle measurable with the 0 analysis and 
not just A 2 and A 3 , as Equation ( 6 ) suggests. In prac¬ 
tice, the A 2 and A 3 are measurable with a greater sta¬ 
tistical significance, because the terms A^/,;(</>) are much 
smaller, for i y 2,3. As a result, these terms affect less 
the 0 distribution. Figures 9 and 10 show the /,; func¬ 
tions for electron acceptances and efficiencies. The shape 
of the fi functions is almost identical for the muons. For 
perfect acceptance and no kinematic cuts (ae = 1 ), the 
only surviving /, functions would be /_ 1 , / 2 , / 3 , f$, and 
fr, and they would be equal to 8/3, 2/3 cos 20, 7 r/ 2 cos</>, 
4/3 sin 20 and 7 r/ 2 sin 20 , in accordance to Equation ( 6 ). 
Figure 11 shows the 0 distribution (9) for the low qx bin, 
with the background neglected and with only one coeffi¬ 
cient varying at a time. We see that the 0 distribution is 
primarily sensitive to A 2 and A 3 , and these coefficients 
are the easier measurable ones with the 0 analysis. Fig¬ 
ure 12 shows the Monte Carlo expected experimental 0 
distributions for a data sample of the size of the Tevatron 
Run I. 




FIG. 11: The expected charged lepton 0 distribution in the 
Collins-Soper W rest-frame for the first qr bin, if we vary 
only one coefficient at the time (with a step of 0 . 1 ) and keep 
the rest at the Standard Model value. Only the A 2 and A 3 
significantly affect the shape of the distributions. The same 
is true for the higher qr bins. 


The final step is to extract the angular coefficients us¬ 
ing the data of Figure 12. We keep the Aj ^ 2 ,3 coefficients 
frozen at their Standard Model values we determined 
above and we fit the distributions to the fi varying A 2 
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FIG. 12: The Monte Carlo prediction for the <f> distribution, 
including the detector acceptance and efficiency, for the four 
qr bins (points). The size of the sample and error-bars cor¬ 
respond to the expected yields at the Tevatron at yfs = 1.8 
TeV). Also shown is the result of the fit (dashed line) using 
the theoretical prediction of Equation (9) and varying A 2 and 
A 3 simultaneously. The backgrounds are ignored, only VF+jet 
events are considered. 

and A 3 simultaneously. The result of the fit can be seen 
in Figure 12 and the extracted coefficients in Figure 13. 
We conclude that the measured angular coefficients are 
close to the values we extracted in section II, verifying 
that the method is self-consistent and could be used for 
an experimental measurement of the W angular coeffi¬ 
cients. The same technique can be applied in Z boson 
experimental studies - which do not demonstrate any 
problems in the kinematic reconstruction of the boson - 
using the future statistically significant datasets of the 
Tevatron and the LHC. 


V. SUMMARY 


FIG. 13: The extracted angular coefficients, using the Monte 
Carlo data of Figure 12 and fitting using the fi functions. 
The angular coefficients that do not affect the shape of the 
4> distribution are kept at their Standard Model values. The 
errorbars are statistical only. 


application of quality cuts to reduce backgrounds. Two 
additional issues are the W mass width effect and the 
resolution of the two-fold ambiguity in the longitudinal 
momentum of the neutrino. We presented the effect of 
these factors on the angular distributions and noted that 
both problems do not affect the azimuthal angle of the 
charged lepton in the CS frame. Finally, we suggested 
a method of extracting the angular coefficients without 
having to divide the experimental data by Monte Carlo 
distributions of isotropic W decays. Passing the genera¬ 
tor data through a detector simulator and analyzing the 
resulting data, we were able to get back the angular co¬ 
efficients we determined from the direct analysis of the 
generator data, demonstrating that this procedure is re¬ 
liable for the experimental measurement of the angular 
coefficients. 


The Standard Model prediction for the angular coeffi¬ 
cients and the associated helicity cross sections of the W 
production in a hadron collider up to order a% in QCD 
and at y/s = 1.8 TeV was presented. The experimental 
measurement of the angular distributions is distorted due 
to the acceptances and efficiencies of the detector and the 
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